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ABSTRACT 



We report on the observation of the bright, long gamma-ray burst, 
GRB 090902B, by the Gamma-ray Burst Monitor (GBM) and Large Area Tele- 
scope (LAT) instruments on-board the Fermi observatory. This was one of the 
brightest GRBs to have been observed by the LAT, which detected several hun- 
dred photons during the prompt phase. With a redshift of z = 1.822, this burst 
is among the most luminous detected by Fermi. Time- resolved spectral analysis 
reveals a significant power-law component in the LAT data that is distinct from 
the usual Band model emission that is seen in the sub-MeV energy range. This 
power-law component appears to extrapolate from the GeV range to the lowest 
energies and is more intense than the Band component both below ~ 50 keV and 
above 100 MeV. The Band component undergoes substantial spectral evolution 
over the entire course of the burst, while the photon index of the power-law com- 
ponent remains constant for most of the prompt phase, then hardens significantly 
towards the end. After the prompt phase, power-law emission persists in the LAT 
data as late as 1 ks post-trigger, with its flux declining as t~^-^. The LAT detected 
a photon with the highest energy so far measured from a GRB, 33.4l3;5 GeV. 
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This event arrived 82 seconds after the GBM trigger and ~50 seconds after the 
prompt phase emission had ended in the GBM band. We discuss the imphcations 
of these results for models of GRB emission and for constraints on models of the 
Extragalactic Background Light. 

Subject headings: gamma rays: bursts 



1. Introduction 



T he Fermi Ganima-r ay Space Telescope hosts two instrun ients, the Large Ar ea Tele- 
scope ( Atwood et al. 2009 ) and the Gamma-ray Burst Monitor ( Meegan et al. 20091 ). which 
together are capable of measuring the spectral parameters of gamma-ray bursts (GRBs) 
across seven decades in energy. Since the start of GBM and LAT science operations in 
early August 2008, emission at energies >100MeV has been detected from ten GRBs. 
These detections were made possible by the LAT's greater sensitivity and shorter dead- 
time (26 /is) compared to previous instruments. Prior to Fermi, high-energy gamma-rays 
from GRBs with energies up to 18 GeV were observed by the EGRET instrument on-board 
the Compton Gamma-ray Observatory. The EGRET observations suggested three types of 
high- energy emission: an extrapolation of the low energy spectra to the >100 MeV band 
(e.g., Dingus et al. 19981), an additional s pectral component during the prompt emission 

7, a GeV after- 



( 1 Gonzalez et al 



20031: 



Kaneko et all |2008|) and in the case of GRB 9402 



glow which was detectable for 90 minutes after the tr igger (IHurlev et al. 



1994f ). The red- 



shifts of these events were not determined. Recently, I Giuliani et al.l (120081 ) reported that 
GRB 0805 14B which triggered AGILE at lower energies, was detected by the GRID instru- 
ment up to 300 M eV. A photometric redshift of 2; = 1.8^*^03 was reported for this event 
faossi et al.ll2008h . 



In the Fermi era, due to the advanced localiz ation capabilities of the LAT and the rapid 
follow-up by the Swift narrow field instruments (IGehrels et al.ll2004l ) and the ground-based 
follow-up community, redshifts for five of the t en LAT bursts have been measured. These 
include GRB 080916C with z = 4.35 ± 0.15 jGreiner et al.ll2009h . a long bu rst that has 



the hi ghest inferred isotropic energy, Eiso ~ 8.8 x 10^^ ergs 
2009dh . and GRB 090510 with z = 0.903 ± 0.003 jRau et al 



( 10 ke V-10 GeV) flAbdo et al. 



20091 ). the second short burst 



seen by the LAT and the first short burst to show defin itively an additiona l hard power-law 
component in the GeV band during the prompt phase (lAbdo et al.ll2009al ). 



GRB 090902B is a long, fairly intense burst with a redshift of z = 1.822 (jCucchiara et al. 



20091 ) and fluence of (4.36 ±0.06) x lO'^ergcm'^ (10 keV-10 GeV) over the first 25 seconds 
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of the prompt emission. These data give an isotropic energy Eiso = (3.63 ±0.05) x 10^^ ergs, 
comparable to that of GRB 080916C. Similar to GRB 090510, GRB 090902B has a significant 
additional, hard power-law component that appears during the prompt phase. Furthermore, 
a spectral feature at energies < 50 keV is evident in the GBM spectrum of GRB 090902B 
that is consistent with an extr apolation of the > 10 MeV power-law emission down to those 



energies. In previous analyses, iPreece et al.l (119961 ) reported evidence for an additional low 



energy spectral component below 20keV for ~ 15% of BATSE bursts. 

We report on the observations and analysis of gamma-ray emission from GRB 090902B 
measured by the GBM and LAT instruments. In section 2, we present details of the detections 
by both instruments and summarize the follow-up observations. In section 3, we show the 
light curves of the prompt emission as seen by the various detectors and describe the extended 
emission found in the LAT data out to 1 ks after the trigger. In section 4, we present the 
time-resolved spectral analysis of the burst emission during the prompt phase. Finally, in 
section 5, we discuss the physical interpretation of the GBM and LAT data, focusing on the 
implications of the power-law component for models of GRB physics. 



2. Observations 



On 2009 September 2 at 11:05:08.31 UT, the Fermi Gamma-ray Burst Monitor trig- 
gered on and localized the b right burst GRB 090202B (trigger 273582310 / 090902462, 
Bissaldi fc Connaughtoru 120091 ). The burst was within the LAT field of view initially at an 
angle of 51° from the boresight. This event was sufficiently bright in the GBM that an Au- 
tonomous Repoint Request was made, and the spacecraft began slewing within 10 seconds 
towards the burst. After ~200 seconds, it had pointed the LAT boresight to within a few 
degrees of the final burst localization. It maintained that pointing until ~1 ks post-trigger, 
when the Earth's limb began to enter the LAT field-of-view (FOV). This burst was detected 
up to ~5 MeV by GBM, and emission was significantly detected by the LAT, with 39 pho- 
tons above 1 GeV. The highest energy photon had E = 33.4^35 GeV and arrived 82 seconds 
after the G BM trigger; and the in itial analyses detected photons as late as 300 seconds after 
the trigger (Ide Palma et al.ll2009l ). 



From the LAT data, the burst was localized to R.A.(J2000), Dec(J2000) = 265.00, 
27.33 with a statistical uncertainty of 0.04° (+ <0.1° systematic), enabling Target of Op- 
portunity observations to begin ~ 12.5 hours after the trigger with the narrow field in- 
struments on Swift. A candidate X- ray afterglow within th e LAT error circle was de- 
tected by the X-Ray Telescope (X RT, Kennea &: Strattalbood ). This source was confirmed 
to be fading (IStratta et al.ll2009l ). and UVOT observations revealed the optical afterglow 
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( ISwenson fc Strattall2009l ). The earlie st ground-based opt ical observations were obtained by 
ROTSE-IIIa ^1.4 hours post tr igger (jPandey et al.ll2009l ). Other detections were repo rted 
in the optical (IPerlev et al.ll2009n. in the near infrared by GRON D (jOlivares et al.ll2009l ) and 
in the radio (Ivan der Horst et al.ll2009l : IChandra fc Frail! 120091 ) . The location of the fading 
source detected by GROND was R A.(J2000), Dec('J200 0) = 17^39'^45!41, +27°19'27.1", 3.3 
arcminutes from t he LAT location (lOlivares et al.ll2009l ). The afterglow redshift of 2; = 1.822 
was measured by ICucchiara et al.l (120091 ) using the GMOS spectrograph mounted on the 
Gemini-North telescope. 



3. Light Curves 



In Figure (H we show the GBM and LAT light curves in several energy bands. The top 
three panels show data from the most brightly illuminated Nal and BGO detectors of the 
GBM, and the bottom three panels show the LAT data with various event selections. In the 
bottom panel, the measured photon energies are plotted as a function of time, including the 
highest energy event [E = 33.4 GeV) that arrived 82 seconds after the GBM trigger time, 
Tq. From the GBM light curves, we see that at energies <lMeV the prompt phase ends 
approximately 25 seconds after Tn. Detailed analysis of the GBM data for energies 50-300 
keV yields a formal T90 duratior|^ of 21.9 seconds starting at To + 2.2 s. By contrast, the 
LAT emission >100MeV clearly continues well after this time range. 

On time scales longer than the prompt phase, the LAT detects emission from GRB 
090902B as late as 1 ks after the GBM trigger. The spectrum of this emission is consistent 
with a power-law with photon index F = —2.1 ± 0.1, and its flux (>100MeV) declines as 
^-i.Sito.i Q^gj, interval (Tq -f- 25, Tq -|- 1000 s). As we note above, the LAT observations 
are interrupted by entry of the Earth's limb into the FOV, but analysis of data after Tq -t- 
3600 s, when the source location is again unocculted, shows that any later emission lies 
below the LAT sensitivity (Figure [2]). The upper limit we obtain for data after Tq + 3600 is 
consistent with an extrapolation of the decay. Similar late time emission for energies 
>100 MeV that extends well beyond the pro mpt phase has been seen for five e arlier bursts by 
Fermz: GRB 0809 16C JAbdo et al.ll2009dh : GRB 090323 Johno et al.ll2009h: GRB 09 0328 
dCutini et al]l2009h: GRB 90510. independently see n bv AGILE dGiuhani et al.l l2009h and 
bv Fermi dohirlanda et"aDl2009h : and GRB 090626 Jpiron et al.ll2009h . 



^The T90 duration is the time over which the central 90% of the counts from the burst have been 
accumulated. 
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4. Time-resolved Spectral Analysis 

Spectral analysis was performed using the data from both the GBM and the LAT. 
These analyses include data from the Nal detectors 0,1,2,9,10 and BGO detectors, and LAT 
"transient" class data, with front- and back-converting events considered separately. The 
Nal data are fit from 8 keV to 1 MeV and the BGO from 250 keV to 40 MeV using the Time 
Tagged Event (TTE) data, which are high time resolution data that allow us to define the 
time intervals based on the structure of the GBM and LAT light curves. The LAT data are 
fit from 100 MeV to 200 GeV. An effective area correction of 0.9 has been fit to the BGO 
data to match the model normalizations given by the Nal data; this correction is consistent 
with the uncertainties in the GBM detector responses. The fits were performed with the 
spectral analysis software package RMFIT (v ersion 3.1). For fu rther details on the dat a 



extraction and spectral analysis procedures see I Ab do et al.l ( l2009bl ) and I Ab do et al.l (l2009d ). 



The time-int egrated spectrum of GRB 090902B is best modeled by a Band function 



( iBand et al.lll993l ) and a power- law component (Tabled]). The power- law component signif- 
icantly improves the fit between 8 keV and 200 GeV both in the time-integrated spectrum 
and in the individual time intervals where there are sufficient statistics. It is also required 
when considering only the GBM data (8 keV-40 MeV) for the time-integrated spectrum, as 
its inclusion causes an improvement of ~ 2000 in the CSTAT statistic over the Band function 
alone. When data below ~ 50 keV are excluded, a power-law component can be neglected 
in the GBM-only fits. We conclude that this power-law component contributes a significant 
part of the emission both at low (<50keV) and high (> 100 MeV) energies. Figure [3] shows 
the counts and unfolded vFi, spectra for a Band function with a power- law component fit to 
the data for interval b (when the low energy excess is most significant) using the parameters 
given in Table [H 

Spectral evolution is apparent in the Band function component from the changing E^^ak 
values throughout the burst, while j3 remains soft until interval e when it hardens signifi- 
cantly. j3 is similarly hard in interval f, after which the Band function component is no longer 
detected. The hardening of (3 is accompanied by an apparent hardening of the power-law 
index, F, which until interval e does not exhibit much variation. However, this is not defini- 
tive since the flux is too low to constrain F in intervals e and f separately. A spectral fit of 
the sum of these two intervals confirms the presence of both a harder (3 and a harder F, with 
a clear statistical preference for the inclusion of the power-law component. An equally good 
fit is obtained in the combined e + f interval if this power-law has an exponential cut-off at 
high energies, with the preferred cut-off energy lying above 2 GeV. Finally, we note that in 
interval b, a marginally better fit is achieved using a model with the additional power-law 
component having an exponential cut-off at high energies. The improvement is at the ~ 3(T 
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level and indicates weak evidence for a cutoff in the second component, placing a lower limit 
on the cutoff energy in this interval of about 1 GeV. 



5. Discussion and Interpretation 



The Fermi data for GRB 090902B show for the first time clear evidence of excess 
emission both at low energies (< 50 keV) and at high energies (>100MeV), while the Band 
function alone fits data at intermediate energies adequately. These excesses are well-fit 
by a single power-law component suggesting a common origin. This power-law component 
accounts for fa 24% of the total fluence in the 10 keV-10 GeV range, and its photon index 
is hard, with a value ~ —1.9 throughout most of the prompt phase. Such a hard component 
producing the observed excess at low energies is difficult to explain in the context of leptonic 
models by the usual synchrotron self-Compton (SSC) mechanisms. 

In the simplest versions of these models, the peak of the SSC emission is expected 
to have a much higher energy than the synchrotron peak at MeV energies, and the SSC 
component has a soft tail that is well below the synchrotron flux at lower energies and 
so would not produce excess emissio n below ~ 50 keV. Hadronic models, either in the 



form of proton syn chrotron radiation ( IRazzaque et al.ll2009l ) or photohadronic interactions 
( lAsano et al.ll2009l ). can produce a hard component with a similar low energy excess via 
direct and cascade radiation (e.g., synchrotron emission by secondary pairs at low energies). 
However, the total energy release in hadronic models would exceed the observed gamma-ray 
energy of Eiso = 3.63 x 10^^ ergs significantly and may pose a challenge for the total energy 
budget. Collimation into a narrow jet may alleviate the energy requirements, since the actual 
energy release from GRB 090902B can be smaller by a jet beaming factor > I/Fq from the 
apparent isotropic value, where Tq is the bulk Lorentz factor of the fireball. 

From the observation of a II.IGIq gg GeV photon in interval c, the highest energy during 
the prompt phase and thus the most constraining, we derive a minimum value of the bulk 
Lorentz factor Fmin ~ 1000 using the flux variability time scale of ty 53 ms found in the 
BGO data. This limit follows from the constraint that the opacity for pair production with 
target photons fit ted by the Band-|-PL model in interval c is less than unity for the 11.1 6 GeV 
photon (see, e.g.. lFenimore et al.lll993l : iBaring fc Hardinpl997l : iLithwick fc Sarill2001l). This 
high F min value is of the same order as the values derived for GRB 080916C ( lAbdo et al. 



2009d|) and GRB 090510 flAbdo et al.ll2009af ). both of which have been detected at > 10 GeV 
with the LAT. 



The delayed onset of the >100 MeV emission from the GBM trigger has been modeled for 
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GRB 080916C as arising from proton synchrotron radiation in the prompt phase (IRazzaque et al 
20091 ) and f or GRB 090510 as arising from electron synchrotron ra diation in the early after- 
glow phase (IKumar fc Barniol Durarul2009l : iGhirlanda et al.ll2009l ). In order to produce the 
peak of the LAT emission at ~ To + 9 s in the early afterglow scenario for GRB 090902B 
from deceleration of the GRB fireball, a value of Tq f« 1000 is required. This is similar to 
Fmin that we calculate, but the observed large amplitude variability on short time scales 
90 ms) in the LAT data, which is usually attributed to prompt emission, argues against 
such models. Also, the appearance of the power-law component extending down to ~ 8 keV 
within only a few seconds of the GRB trigger disfavors an afterglow interpretation. The 
proton synchrotron model, on the other hand, requires a rather large total energy budget, 
as mentioned previously. 

Yet another interpretation of the observed excess in the high and low energies may 
be provided by two non-t hermal power-law comp o nents along with a thermal component 
from the jet photosphere (IMeszaros fc Reed l2000l : iRydd |2004| ). The thermal component, 
broadened by t emperature varia tions, then accounts for the > 100 keV-few MeV emission 
with To ~ 930 (jPe'er et al.l 120071 ). although fits of such a model to our data do not improve 
over the Band-f-PL model. Furthermore, it is difficult for the photospheric model to explain 
the delayed onset of the > 100 MeV emission. 

The detection of the 33.4 GeV photon, 82 seconds after the GRB trigger and well after 
the soft gamma-ray emission subsided, may help constrain the origin of the late-time decay 
of the power-law component, which goes as A synchrotron origin of the 33.4 GeV 

photon would be difficult since it would require significant energy gain by electrons over a 
gyroradius and a bulk Lorentz factor > 1500. In the case of diffusive shock-a cceleration, the 
energ y losses in the upstream region of the shock may dominate (see, e.g., iLi fc Waxman 
20061 ) and prevent acceleration of electrons to an energy high enough to radiate a 33.4 GeV 
photon. An interpretation by afterglow SSC emission is still possible, however. 

The cons traints on the quantum gravit y mass scale from GRB 090902B using the time- 
of-flight test ( Amelino-Camelia et al.lll998l ) are much weaker than those from GRB 090510 



( lAbdo et al.ll2009a( ) due to the larger interval, 82 seconds, between Tq and the arrival time of 
the 33.4 GeV photon. However, the moderately high redshift {z = 1.822) of GRB 090902B 
allows us to use this photon to probe and constrain models of the Extragalactic Background 



2008; Gilmore et al 



2009 



Light /EBL: iKneiske et al.ll2004l: Metcalfe et al.ll2003l : IStecker et al.ll2006l : IPranceschini et al. 



Pinke et al. 



20091 ). The 33.4 GeV photon would not be absorbed 



by the EBL i n any models except for the "fast evolution" and the "baseline" models by 
Stecker et al.l (120061 ) . which give optical depths of r^^ = 7.7 and 5.8, respectively. We have 
performed spectral fits of the LAT data with and without the predicted EBL absorption 
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from Stecker's models assuming a simple power-law as the intrinsic emission model. Based 
on Monte-Carlo simulations, we found that Stecker's fast evolution and baseline models are 
disfavored at a > 3cr level. 

In summary, GRB 090902B is one of the brightest bursts detected by the GBM and 
LAT instruments on Fermi. It clearly shows excess emission at high and low energies during 
the prompt phase, requiring a hard power-law component in addition to the usual Band 
function in order to fit the data. The origin of this component is not understood, and its 
presence in this burst poses genuine challenges for the theoretical models. Like the other two 
bright Fermi bursts detected by the LAT, GRB 080916C and GRB 090510, GRB 090902B 
appears to possess a very high Lorentz factor for the bulk outflow, F ^ 1000, and has some 
suggestion of a delayed onset of the emission above ~ 100 MeV. Finally, the 33.4 GeV photon, 
the highest energy yet detected from a GRB, and the z = 1.822 redshift of this burst have 
allowed us to place significant constraints on some models of the Extragalactic Background 
Light. 
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Fig. 1. — GBM and LAT light curves for the gamma-ray emission of GRB 090902B. The 
data from the GBM Nal detectors were divided into soft (8-14.3 keV) and hard (14.3-260 
keV) bands in order to reveal any obvious similarities between the light curve at the lowest 
energies and that of the LAT data. The fourth panel shows all LAT events that pass the 
on-board gamma filter, while the fifth and sixth panels show data for the "transient" class 
event selection for energies > 100 MeV and > 1 GeV, respectively. The vertical lines indicate 
the boundaries of the intervals used for the time-resolved spectral analysis. Those time 
boundaries are at Tq + (0,4.6,9.6,13.0,19.2,22.7,25.0,30.0) seconds. The insets show the 
counts for the corresponding dataset binned using these intervals in order to illustrate the 
relative numbers of counts considered in each spectral fit. 
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Fig. 2.— Light curve of GRB 090902B for energies 0.1-300 GeV from unbinned likelihood 
fits to the LAT data. After the prompt phase, extended or afterglow emission consistent with 
a temporal profile oc t^^'^ (dashed line) lasts until ~ Tq + 1000 s. The upper limit at times 
> To + 3600 s was derived from the data collected after the source emerged from occultation 
by the Earth. 
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Fig. 3. — Joint fit of GBM and LAT data to interval b, (Tq + 4.6, Tq + 9.6 s). Top: Counts 
spectrum; separate model components are plotted, Band (dashed), power-law (solid). Bot- 
tom: Unfolded vF^ spectrum. The extension of the > 100 MeV power-law component to the 
lowest energies (< 50keV) is shown. 



Table 1. Band function + power-law fit parameters for the time-resolved spectral fits. 



Interval 


Time Range (s) 


Epeak (kev) 


a 


13 


r 


CSTAT/DOF 


ACSTAT 


Energy fluence 
(erg cm-2, 8 keV-30 GeV) 




0.0-30.0 


726 (±8) 


-0.61 (±0.01) 


-3.8 (1^:^) 


-1.93(tSri) 


2562/963 


2005 


(4.59 ± 0.05)xl0~* 


a. 


0.0-4.6 


526 (±12) 


-0.09 (±0.04) 




-i.87(«;°t) 


901/963 


43 


(3.72 ± 0.13)xl0-5 


b. 


4.6- 9.6 


908 (til) 


0.07 (±0.03) 


-3.9 n.i) 


-1.94 (±0.02) 


1250/963 


3165 


(1.44 ± 0.03)xl0-* 


c. 


9.6-13.0 


821 (±16) 


-0.26 (±0.03) 




-1.98 (±0.02) 


1310/963 


2109 


(9.42 ± 0.24) xlO"'' 


d. 


13.0-19.2 


529 (±9) 


-0.65 (±-0.02) 


-3.2 ni) 


-1.86 (± 0.02) 


1418/963 


199 


(1.29 ±0.03) X 10-* 


e. 


19.2-22.7 


317 (±8) 


-0.78 (±-0.02) 


-2.4 (±0.1) 




1117/965 




(4.8 ± 0.2) xlO-'^ 


f. 


22.7-25.0 


236 (til) 


-1-30 W 


-2.2 (±0.1) 




1077/965 




(1.0 ± 0.1)xl0-5 


e.+f. 


19.2-25.0 


327 (±8) 


-0.91 (±0.02) 


-2.6 (±0.1) 


-1.59 (±0.20) 


1219/963 


16 


(6.1 ±0.4)xl0-5 


g- 


25.0-30.0 








-1-93 if^il) 


1209/967 




(6.8 ± 0.8)xl0-® 



Note. — The time range values are relative to the trigger time Tq. The column ACSTAT gives the change in CSTAT when fitting with only the Band 
function versus Band+power-law. The Band function is given by 

^ ' VlOOkeVy "^V iJpoak J 

' exp(/3 I , E > Ec, (1) 



100keV(2 + a)y \ lOOkeV 



where Ec = (a - /3)£;pcak/(2 + a) llBand et alj|l993h . The power-law function is given by 



